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We investigated a nuclear modification difference between up- and down- valence quark distributions by analyz- 
ing structure function F2 and Drell-Yan cross-section ratios. Although nuclear modifications of the valence-quark 
distributions themselves are rather well determined, it is difficult to find their difference from the present data. 
We estimated such an effect on the NuTeV sin'^ 9w value and its uncertainty by the Hessian method. At this 
stage, it is not large enough to explain the whole NuTeV anomaly. However, the modification difference cannot 
be precisely determined, so that further studies are needed. 



1. Introduction 

Weak-mixing angle sin^ 9w is one of the funda- 
mental constants in the standard model. It has 
been measured experimentally by various meth- 
ods such as atomic parity violation and polar- 
ization asymmetry in electron-positron annihi- 
lation. The NuTeV collaboration determined 
it by neutrino- and antineutrino-nucleus scat- 
tering and found that its value, sin^ 6w = 
0.2277 ± 0.0013 (stat) ± 0.0009 (syst) in the on- 
shell scheme, is significantly larger than other 
measurements (sin^ = 0.2227 ± 0.0004) [T]. 
It is called "NuTeV sin^ 6w anomaly" . 

Because there may be new physics behind this 
anomaly, careful analyses are needed for clarify- 
ing the situation. In our work, we investigate a 
conservative explanation without exotic mecha- 
nisms. The NuTeV target is iron, so that nuclear 
corrections could be a candidate for the anoma- 
lous result. There are various factors including 
the effects of neutron excess, strange-antistrange 
asymmetry, isospin violation, and modification of 
valence-quark distributions. Among them, the 
nuclear modification difference between Uv and 
dy is discussed in this paper. It was first investi- 
gated in Ref. 2 and a detailed analysis has been 
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done in Ref. [H]. Here, we report recent analysis 
results. 

This paper consists of the following. In section 
121 our analysis method and results are explained 
for determining the nuclear modification differ- 
ence. Its effect on sin^ 9w is calculated in section 
131 The results are summarized in sectional 

2. Nuclear modification of valence-quark 
distributions 

It is known that nuclear parton distribution 
functions (PDFs) are modified from those of the 
nucleonic distributions. Such nuclear modifica- 
tions have been investigated especially in the 
structure function F2 . Determination of each par- 
ton modification is not straightforward from F2 
and Drell-Yan data. However, gross features of 
nuclear PDFs are now determined, for example, 
in Ref. 4 . Among the nuclear PDF corrections, 
it was pointed out that valence-quark modifica- 
tions affect the sin^ 9w determination [2| . 

In order to discuss such a nuclear effect, we 
define nuclear modification factors Wu^ and w^^, 
for up- and down- valence quark distributions by 

^ Zuy{x) +Ndy{x) 
(x) = Wu, (x. A, Z) -J , 

,Ai \ t A 7\ Zdy(x) + Nuy{x) 
d^{x) = Wd^ [x, A, Z) , (1) 
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where and d^{x) are the distributions in 

the nucleon, Z is the atomic number of a nu- 
cleus, and A is the mass number. For simphcity, 
the dependence is abbreviated. It should be 
noticed that these relations are used at any 
for the present research, although similar equa- 
tions are defined only at Q^=l GeV^ in Ref. 0]. 
These equations are motived by the following con- 
siderations. If there were no nuclear correction, 
the distribution of a nucleus is given by the 
simple summation of proton and neutron contri- 
butions Zu^ + Nu". The isospin symmetry is 
usually assumed for the parton distributions, so 
that it becomes Zuy + Ndy. It is divided by A 
{{Zuy + Ndy) /A) because we use nuclear PDFs 
per nucleon. Therefore, the function it;„^, indi- 
cates nuclear correction to this distribution. 

In particular, we are interested in the modi- 
fication difference Wu^ — Wd^_ and its effect on 
the sin^ 9w determination. Therefore, specific 
parameters (a^, b'^, c^, d^) are assigned for the 
difference at Q^^l GeV^ 0: 
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Because of the baryon-number and charge conser- 
vations, the parameter a'^ is fixed, and then the 
number of parameters is three. The parameter 
/?„ and other parameters in Wu^ + Wd„, wq, and 
Wg are fixed at the values of our recent analysis 
[3]. The mass- number dependence is assumed in 
the 1 — 1/A^/^ form simply by considering nu- 
clear volume and surface contributions The 
X dependence is motivated by the shape of nu- 
clear modification data of F2. However, we should 
aware that appropriate A and x dependence is not 
known almost at all for Wu^ — Wd^ ■ 

We determined the parameters by using exper- 
imental data for the nuclear F2 ratios and Drell- 
Yan cross-section ratios as investigated in Ref. 

The parameters 6^, c^, and d'^ are optimized 
by a analysis with the data. The obtained 
distribution is shown by the ratio 
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Wd^ {x, A, Z) - Wu„ [x, A, Z) 
Wd^ [x, A, Z) + Wu„ {x, A, Z) 
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Figure 1. The obtained distribution ey{x) is 
shown at Q^=20 GeV^ [HI with previous results 
P]. The details are explained in the text. 

in Fig. n Equation Q is used at any Q^. The 
distribution ey{x) is shown by the solid curve, and 
the shaded area indicates the one-cr error range. 
The error is estimated by the Hessian method by 
using the determined parameters and error ma- 
trix in the analysis. The distribution is shown 
at (52=20 GeV^ which is about the average (J 
value of the NuTeV data. 

The dashed and dotted curves indicate the esti- 
mated distributions in Ref. 2 . The dashed curve 
is obtained by calculating — —{N — Z){uy — 
dv){'Wy — l)/[A{uy + dy)wy], which is one of the 
candidates for satisfying the baryon-number and 
charge conservations. The function is given 
by Wy — {wu^ -\- Wd^)/2. The dotted is obtained 
by Eq. with the 2001 version of the nuclear 
PDFs , 4 . As obvious from the figure, three dis- 
tributions are much different. However, they are 
well within the error band, which indicates that 
these results are consistent each other. 



3. EfTects on NuTeV sin^ 
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From the neutrino- and antineutrino-nucleon 
scattering, sin^ 9w could be obtained by using the 
Paschos-Wolfenstein (PW) relation: 
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where and a 



are charged-current (CC) 
and neutral-current (NC) cross sections. There 
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are various nuclear corrections to the PW rela- 
tion. Expanding the expression in terms of the 
correction factors, we obtain a modified PW re- 
lation for a nucleus [2]: 

Ra = ^- sin^ 9w 

- £.(a;)| Q - sin^ Ow^ | ^ j| _ j^jz " I ^in' 0w| 
+ 0(e2) + 0(e„) + 0(e,) + 0(eJ. (5) 

The first correction term with £i,(x) is investi- 
gated in this paper. In addition, there are cor- 
rections associated with the neutron-excess (£«), 
strange-antistrange asymmetry {Ss), and charm- 
anticharm asymmetry (e^) factors. The variable 
y is defined by y — q'~' / with the energy transfer 

Because the e.u correction term depends on the 
variables a:, y, and Q^, we need to take the NuTeV 
kinematics into account to estimate an effect on 
the sin^ 6w determination. In particular, there 
are few experimental data in the large x region, 
where the distribution ey{x) becomes significant 
as shown in Fig. ^ Namely, the Sy effect on 
sin^ 6w could be much smaller than the one ex- 
pected from Fig. ^ Fortunately, such kinemat- 
ical factors are supplied in Ref. 5 . Comparing 
our PDF definition with the NuTeV one, we find 
the relations 

6ul = u*p - d*„ = -Sv (w„„ -I- Wd^)xuy, 

6dl = dip - M*„ = +ey {wu^ + Wdjx dy, (6) 

where the asterisk (*) indicates the NuTeV dis- 
tributions. Therefore, the nuclear modification 
difference £„ corresponds to the isospin violation 
in the NuTeV terminology. Then, the contribu- 
tion to the NuTeV sin^ 9yir is given by 

A(sin2 ew) = ~ J dx { F[5ul,x\ 5ul{x) 

+ F[5d:,x]5dl{x)], (7) 

where F[6u*,x\ and F[5d*^x\ are the functionals 
provided in Fig. 1 of Ref. [HI. Our sign con- 
vention of A(sin^6'vi/) is opposite to the NuTeV 
one in Eq. Q. The distributions 5ul{x) and 



5d*{x) are calculated by Eq. © with the distri- 
bution Ey {x) , which was obtained in the previous 
section. Calculating the integral numerically, we 
obtain 3 

A(sin2 Ow) = 0.0004 ± 0.0015, (8) 

at Q2^20 GeV2. 

In comparison with the NuTeV deviation 

0. 0050. the correction is an order of magnitude 
smaller. However, the error becomes comparable 
to the deviation. The magnitude of the error de- 
pends much on the analysis condition. Therefore, 
careful analyses could alter the values in Eq. . 
In particular, the nuclear modification difference 
Wu^ — Wd^ , which we have investigated in this pa- 
per, is not a uniquely determined quantity at this 
stage. In order to determine this difference, we 
need future experimental efforts. 

4. Summary 

We have analyzed the experimental data on nu- 
clear structure functions F^ and Drell-Yan cross 
sections for extracting the nuclear modification 
difference between Uy and dy distributions. Us- 
ing the obtained distribution, we investigated its 
effect on the NuTeV sin^ 9\y determination. We 
found a rather small effect; however, the uncer- 
tainty on the sin^ 6iy is not small in comparison 
with the NuTeV deviation. Because such as nu- 
clear modification is difficult to be determined at 
this stage, we need further studies to pin down 
the modification. 
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